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For the case of phase damping (pure decoherence) we investigate the extent to which environmental 
traits are imprinted on an open quantum system. The dynamics is described using the quantum 
channel approach. We study what the knowledge of the channel may reveal about the nature of its 
underlying dynamics and, conversely, what the dynamics tells us about how to consistently model 
the environment. We find that for a Markov phase-damping channel, that is, a channel compatible 
with a time-continuous Markovian evolution, the environment may adequately be represented by a 
mixture of only a few coherent states. For arbitrary Hilbert space dimension A'^ > 4 we refine the idea 
of quantum phase damping, of which we show a means of identification. Symmetry considerations 
are used to identify decoherence-free subspaces of the system. 

PACS numbers: 03.65.Yz,02.50.Ga,03.65.Aa 



I. INTRODUCTION 

Decoherence describes the loss of characteristic traits 
of quantum theory. For the success of emerging quantum 
technologies a detailed understanding of decoherence is 
of great relevance. Schemes to avoid and counter its ef- 
fects need to be developed. Besides, decoherence offers 
insight into the much-debated quantum-to-classical tran- 
sition [ll-Il]- The microscopic dynamics leading to deco- 
herence might be based on very diverse grounds, reaching 
from purely classical phase kicks to a quantum mechani- 
cal formulation based on coupling the system of interest 
to some quantum environment. Hence, a further charac- 
terization of different microscopic mechanisms leading to 
decoherence is desirable. 

Phase damping (or dephasing) denotes the case of pure 
decoherence, corroding the coherences of a quantum state 
while leaving the probabilities, that is, the diagonal ele- 
ments of the density matrix, intact. The dissipation-less 
transition of a pure state into a classical mixture when 
described in the basis of energy eigenstates may serve as 
an example. Despite its simple nature, phase damping is 
enough to completely disentangle quantum states [4] . 

For weak system-environment coupling and short en- 
vironmental correlation times decoherence may be mod- 
elled in terms of Markovian dynamics Q. Here, the fu- 
ture evolution depends solely on the system's present 
state, rather than on anterior times. Yet, there are of 
course instances where this approximation is not valid. 
Given the dynamics of a quantum system it would be 
valuable to have a means of deciding whether the dy- 
namics is Markovian or not. This point has been studied 
lately both in a continuous approach based on the infor- 
mation flow between system and environment Q , as well 
as from a snapshot point of view 0, [1] , where the contin- 



uous dynamics is by construction unavailable. Rather, 
the state of the quantum system is known at separate 
times, only. 

Another interesting question is whether the phase 
damping is due to coupling to a "real" quantum- 
mechanical environment, or wether it can equally be ex- 
plained in terms of stochastically fluctuating, classical 
fields • The latter is a convex combination of uni- 

tary transformations, that is, random unitary (RU) dy- 
namics. While phase damping of a single qubit or qutrit 
may always be described as RU dynamics, in Hilbert 
spaces of dimension A'^ > 4 one cannot always find such 
a representation 

In the article at hand we study the characteristics of 
phase damping from an environmental point of view. 
Phase damping is described utilizing the overlap of dy- 
namical vectors relative to the phase damping basis. The 
nature of the dynamics is refiected by the set of dynam- 
ical vectors, or, conversely, the properties of the dynam- 
ical vectors determine the dynamics to a certain extent. 
In this context, we show that in case of Markovian phase 
damping the dynamical vectors can be identified with co- 
herent states. Likewise, we give instructions for a phys- 
ical model of "quantum phase damping" for arbitrary 
Hilbert space dimension > 4, that is, phase damping 
which does not allow for a RU representation. 

The article is structured as follows. Section [H] 
overviews the theoretical background and serves as an 
introduction to the formal notation. In Sec. Illll we exem- 
plary study phase damping on a single qubit, where all 
characteristics introduced so far actually coincide. Sec- 
tions IIVI and |V] address the Markovianity and the pos- 
sibility of finding a RU representation, respectively. In 
Sec. |VT] we discuss the appearence of decoherence-free 
subspaces due to symmetries in our formalism. 
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II. QUANTUM CHANNELS 



RU Channels 



Based on the fundamental assumption of no initial cor- 
relations between the system g and its environment, the 
most general quantum evolution is given by a completely 
positive map £ : g t-^ £{g). In a Hilbert space of dimen- 
sion A'', these maps (or "quantum channels" ) can always 
be written in terms of at most N'^ Kraus operators Ki 
such that 



g^ g' ^£{g) = Y,K,gKl 



(1) 



(here and in the following we denote the initial state by 
g and its map by g'). It is usually assumed that the 
map is trace-preserving, '^iKIKi = 1, so as to preserve 
probability. If, in addition, the completely mixed state 
is mapped onto itself: KiXj = 1, the channel is said 
to be unital or doubly stochastic [IJ]. Throughout the 
article we will assume that g and g' live in the same 
Hilbert space, that is, the channel £ maps the set of 
states on a Hilbert space of dimension N onto itself. 

When considering a quantum channel of form ([T]), no 
particular assumptions are made about the nature of the 
underlying continuous dynamics. Rather, only a snap- 
shot of the quantum system at a given time is revealed. 
Nevertheless, in some cases it is possible to gather in- 
formation about the nature of the physical processes in- 
volved. In the remainder of this section, we want to dis- 
cuss how certain additional assumptions about the struc- 
ture of the channel may set restrictions to the underlying 
dynamics or vice versa. 



Markovian Channels 

A quantum channel £ is said to be Markovian, if there 
exists a generator >C of a quantum dynamical semigroup 
and a time t > such that 



£{g) 



e^'g 



(2) 



0,01 ■ That is, the channel may be understood as a snap- 
shot of a time-continuous Markovian evolution. The gen- 
erator C may be written in Lindblad form [5)] 



1 ^ 

i=l 



The Markov property of a channel is closely related to 
the notion of infinite divisibility 0, [13 • ^ channel £ is 
called infinitely divisible if, for all ly G V, there exists a 
channel £„ with {£,y)'^ — ^- Surely, a Markov channel £ is 
infinitely divisible: for any given g M it can be written 
as i^-fold concatenation of the channels £i, = e'-'^l" . The 
converse statement, however, is not true in general [l5| . 



One of the standard approaches to the quantum chan- 
nel formalism is based on the reduced dynamics of a 
system interacting with its environment [lO|. In this 
context, decoherence of an open quantum system is in- 
evitably linked to growing entanglement between system 
and environment [2[. Yet, there are instances of irre- 
versible dynamics that may be modeled entirely without 
invoking a quantum environment at all. An important 
example is given by RU dynamics, where the quantum 
channel may be written as a convex combination of uni- 
tary transformations 



> 0,^K 
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The dynamics may thus be thought of as originating from 
classical fluctuations, hence also termed "random exter- 
nal fields" [13] . It is known, for example, that for a 
sing le qubit all doubly stochastic channels are of RU type 
These RU channels gain some significance in the 
field of quantum error correction, where they stand out 
due to the fact that they may be undone completely [3 ■ 
More recently they have also been applied to quantum 
networks [T7| . 



Phase-Damping Channels 

Phase-damping channels are among the simplest con- 
ceivable quantum channels. They are defined by the 
requirement that in a given basis {\n)} — the phase- 
damping basis — no population transfer takes place. The 
only effect of the "environment" is thus to change coher- 
ences {n\g\m) with n ^ m and to leave all {n\g\n) with 
n = I, . . . ,N untouched. In other words, the projectors 
are constants of motion: £(\n){n\) — \n){n\ for all n. 

We conclude that the Kraus operators have to be di- 
agonal in this basis, Ki =diag(aii, 0^2, . • . , aiw) and, cor- 
respondingly, the whole map £ is diagonal, too. We find 



(4) 



with {a„ = (ai„, a2m • • ■ , a™)} any set of N normalized 
complex vectors. It is then sometimes convenient to in- 
troduce the matrix D with Dmn — {o-n\o-m) to write the 
phase-damping channel in the short form g' = D * g, 
where * is the Hadamard product^ that is, the entry-wise 
product of matrices of the same size: g'^„ = DmnQmn 
[18| . From these considerations it is clear that phase- 
damping channels are among the doubly stochastic chan- 
nels. 

If the quantum channel £ is defined via the system's 
coupling to a quantum mechanical environment, the vec- 
tors \an) may be seen as relative states of the environ- 
ment, that is, relative to the states of the distinguished 
basis (see also Sec. |V]). Then the overlap (a„|am), seen 
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as a function of time, may be related to studies of fi- 
delity decay [l^. Yet, this relative state picture need 
not hold in general: the case of RU dynamics shows that 
in certain circumstances decoherence may be attributed 
to stochastic, fluctuating "classical" fields. 

In many situations the dynamical vectors |a„) are of 
course unknown a priori. In particular this holds true in 
an experimental setup where the matrix D is acquired via 
quantum process tomography [lol | . One way of obtaining 
dynamical vectors |a„) from D is by using the Cholesky 
factorization (20| . Given the non- negative matrix D, the 
Cholesky factorization gives D = LL\ with L a lower 
triangular matrix (L is in general not unique). The n-th 
row of L may then be identified with a complex vector 
|a„) € such that = {0'm\0'n) ■ If -D is a positive 

semi-definite matrix of rank r < d, there exists a unique 
L with columns r + 1 through d identical to zero [20j. 
That is, the vectors |a,i) may be chosen as elements of 
C. In the following sections we will study what these 
dynamical vectors \an) reveal about the nature of the 
underlying dynamics. 



III. THE SINGLE QUBIT CASE 

Without revealing too much about the details we want 
to state some results of the following sections. The case 
of a single qubit stands out due to the fact that a phase- 
damping channel is always Markovian (i.e., a snapshot 
of Markovian dynamics) and it is of RU type. These 
findings of course do not allow for generalization to higher 
dimensional systems, yet they have some potential for 
building intuition. For a more rigorous approach as well 
as some missing definitions see Sees. IIVI and IVl 

For a single qubit the phase-damping map is defined 
by the matrix 

D=( , ] , ^^t^M. (5) 

Thus, a single complex number (a2|ai) —: c with modu- 
lus less than one determines the most general single-qubit 
phase- damping channel. Infinite divisibility of a phase- 
damping channel has to be formulated in terms of the 
Hadamard product (see also Sec. IIV|) . that is, the ma- 
trix with {Dy)ran = {DranY^^ has to be checked as 
to its positivity. It is quite straightforward to see that 
the matrix D in Eq. ^ passes this test, which lets us 
conclude that a single qubit phase-damping channel is 
always Markovian (see also Sec. lIVp . 

Another remarkable feature of single qubit phase 
damping — which we will later show to be intimately 
connected to Markovianity — is that the dynamical vec- 
tors in Q may be chosen from the set of coherent 
states {|a) \a G C} of a harmonic oscillator. These are 
eigenstates of the annihilation operator, a\a) — aja), 
and may be seen as displaced vacuum states: \a) ~ 

^aa^-a'a |q^ ^ j^^ ^^^^^ ggg ^^lis, UOtC that for C 7^ 



we may simply let 

c =: e-^^-"^ (6) 

with 7 e R-|_ and oj E [0, 27r). We then define the two- 
mode coherent states |a„) := e^'"" |V7^n)i where h — 
(1,1), ^2 = (1,-1) and cji — ~UJ2 — cj/2. It is easy to 
see that these states give the correct overlap, that is, 
(q!2|q!i) = c. In this vein we can thus always define the 
channel in terms of coherent states |q;„) , n = 1, 2, leading 
to Markovian dynamcis. Written in Lindblad form the 
master equation attains its well-known form 

= -i^[crz,g]-'^[crz,Wz,g]], (7) 

where channel ([5]) with c from ^ is obtained as a snap- 
shot for t — 1. 

Alternatively, we may choose to write the overlap of 
states in the form (a2|ai) = (2p— l)e^'^ with < p < 1, 
obtaining the common quantum channel representation 

m 

E{q) = e-'^'^^ {pQ + (1 - p)a,_Ba,) e^^^^ . (8) 

In this notation it is rather obvious, hence, that the chan- 
nel is RU, which is always true for a single qubit or qutrit 

mm- 

IV. PHASE-DAMPING MARKOV PROCESSES 
AND COHERENT STATES 

As introduced in Sec. |lll the mapping of a quantum 
state subject to phase damping may be written using the 
Hadamard product. Infinite divisibility is equivalent to 
positivity of the matrices D^^ where {D^Y — * . . . ★ 
= u € i.e., {D,)^n = {D^nY^"- While it 
is clear that every Markov channel is infinitely divisible, 
note that the converse also holds in case of phase damping 
(when aU Dmn 7^ 0). 

The argument is based on a theorem by Denisov [l^ 
which states that an infinitely divisible channel £ is of 
the form £ — e^E, where E is an idempotence with 
ECE = CE 0, Ell. In case of phase damping, however, 
the diagonal character of the map together with the re- 
lation E£ = Ee^E = e^E = £ implies E = 1 already, 
whenever all Dmn 0. Thus, Markovianity follows di- 
rectly from infinite divisibility in this case. Recall that 
this is certainly not true for channels in general. 

From Sec. IIIII we already know that any single-qubit 
phase-damping channel (with c ^ 0) is infinitely divisi- 
ble and hence Markovian, but what about higher dimen- 
sions? A simple example shows that already for a 3- 
dimensional quantum system positivity may be violated: 
Let a 3-state phase-damping channel be given by 

(1 ia —ia \ 
—ia 1 a , 
ia a 1 I 
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with real a. Then for ^ < a < 5 the matrix D is posi- 
tive, but all (Hadamard) square roots D^^^ have one neg- 
ative eigenvalue equal to 1 — %/3a. We have thus found 
a single-qutrit phase-damping channel which may not be 
identified as a snapshot of Markovian evolution. 

The notion of infinite divisibility in this section im- 
plicitly assumes all fractional powers of the initial chan- 
nel to be phase damping and therefore diagonal. This 
excludes the rather peculiar case where some dynamics 
is phase damping for a particular time t only, but may 
well change populations at other times. Consider, for 
example, the unitary one-qubit map £{q) = UgU^ with 
U ~ cxp(— zTTcr^t), which is trivially phase damping for 
t= 1,2,3,.... 

From these considerations it is clear that the under- 
lying generator £ of the Markov dynamics is diagonal: 
C{\m){n\) — Zmn\iTi){n\. This, in turn, assures the diago- 
nal character of both the Hamiltonian H =: '^n\n){n\ 

and the Lindblad operators Li =: X]n^"^l")("l (^^^ 
App. thereby leading to the relation 

Zmn = -i{i^m " W„) + {ln\lrn) - ^(||'m|P + ll^nlT), (9) 

where := (ln \ ■ ■ ■ , In"^), r is the number of Lindblad 
operators in ([3]). 

An M-mode coherent state may be written as a dis- 
placement of the vacuum 



M 



with bosonic creation operators [21j . For two coherent 
states I a) this leads to an overlap 



_ AP,a)-^{\\af + \\l3f) 



(10) 



where (•, •) denotes the standard scalar product in C^^ 
and ||a|p = {a, a). 

Based on comparison of Eqs. © and ([TU| we define 
r-mode coherent states of the following form 



-iuJn t 



In Vi 



and find that for Markovian phase damping we may de- 
fine the channel in terms of coherent states such that 



\t=l — x^nvv It^r 

Any Markovian phase-damping channel £ may therefore 
be obtained as the reduced dynamics of the system inter- 
acting with an environment of harmonic oscillators, all in 
coherent states. At first sight, the time dependence of the 
coherent states |a„(t)) = g-'^^rit jy^;^^ niay seem quite 
queer. Yet, this should not be too surprising given that 
a finite reservoir would normally lead to memory effects. 
In order to preserve Markovianity the dynamics thus has 
to be strongly driven, so as to prevent the back-flow of 
information from the environment to the system Q . 

As a final remark we add that the v^-dependence of 
the centroid of the environmental coherent states refiects 
the fundamental relevance of standard Brownian motion 
for all (continuous) Markov processes. 




FIG. 1: (Color online) The construction of an extremal phase- 
damping channel is based on a bipartite system of qudits A 
and B, locally coupling (via k"^ and k^) to a qudit reservoir 
7^. 



V. RU VS. QUANTUM PHASE DAMPING 

We have seen in the previous section that for qutrits — 
or larger systems — phase-damping channels need not be 
Markovian. In a similar spirit, one may ask whether RU 
representations exists for any dimension: can all phase- 
damping processes be written as a convex sum of unitary 
maps? For a qubit, as for the question of Markovianity, 
the answer is positive. In general, however, the answer is 
no, as can be found in fll|. For a two-qubit system, that 
is TV = 4, a physical model for such a non-RU {or quan- 
tum) phase-damping channel is described in Ref. [l3|. 

Our aim here is to offer a method how to construct 
non-RU phase-damping channels in arbitrary dimension, 
extending earlier work. These serve as specific examples; 
it is an entirely different and challenging matter to test 
a given phase-damping channel for the RU property. 

Our method of identification of such a quantum phase- 
damping channel rests on extremality with respect to the 
convex set of doubly stochastic channels. Due to a re- 
sult by Landau and Streater it is known that there exist 
non-unitary, extremal maps in the convex set of diago- 
nal doubly stochastic maps 11']. Extremality is guaran- 
teed for channels where the projectors |a„)(a„| obtained 
from the dynamical vectors {|ai),...,|ajv)} C C in Eq. 
(HI form a (possibly overcomplete) operator basis on C. 
Note that extremality requires r"^ < N [remember that 
r denotes the number of operators used in Eq. ([T} or, 
likewise, the dimensionality of the vectors |a„), N is the 
dimension of the quantum system] . 

The construction of the channel rests on a Hamiltonian 
H locally coupling two qudits (d-dimensional quantum 
systems) to a single qudit environment (cf. Fig.[T]). Then, 
by construction, r'^ = cP = N . In the usual notation we 
set 



H — Hs + Hi + H-n, 



(11) 



where Hs and H-^i denote the Hamiltonian describing 
system and reservoir, respectively. The local coupling of 
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qudits A and B to the reservoir TZ may be set to 



i-.j 



In order to invoke a phase-damping channel on the sys- 
tem we have to require Hg as weU as ah operators ti^, erf 
to be diagonal (the cr-operators will be specified below). 

For any given time t and assuming the usual product 
initial state, g ® a, this dynamics leads to the phase- 
damping channel 

- tm{Uig®a)U^) . (12) 

Due to the restriction to diagonal system Hamilto- 
nian and diagonal coupling, the unitary map U allows 
for a diagonalization in the phase-damping basis {|n)}. 
The interaction may thus be expressed in fashion of a 
controlled-unitary operation 



l/{d^ - 1)! det[ (62-61) (63-61) ••• {bd^--h)] 
spanned by the real vectors 61 , . . . , 6(^2 . In this geometric 
picture we can infer that the channel is extremal iff the 
Bloch vectors 6„ do not point to the same hyperplane in 
R'' or, equivalenty, iff the — 1 dimensional volume 
V spanned by the Bloch vectors is different from zero. 

While still not a general test for the RU property, we 
would like to note that, nonetheless, criterion (jisp may be 
used to give a constructive test of a channel's extremal- 
ity. Given an arbitrary phase-damping channel D, the 
Cholesky factorization gives, as introduced in Sec. [TTl a 
set of dynamical vectors |a„) G C. Recall that r denotes 
the rank of the matrix D. Any r^-dimensional subset of 
the corresponding Bloch vectors 5„ has now to be checked 
for linear independence. If linear independence is found 
in any subset, then — following the equivalence in (1151) — 
we may conclude upon extremality of the channel. For 
r 7^ 1 this immediately excludes the RU property. 



U 



E 

n=l 



n){n\ <Si U„ 



(13) 



Assuming the initial state of the reservoir to be pure, 
that is, cr = |?/'o)(V-'o|: "we obtain the phase-damping chan- 
nel 



= {'^n\i',n) gy, 



(14) 



in terms of the dynamical vectors \tpn 
l,...,d\ 

The properties of the phase-damping channel are now 
encoded in these relative environment states In 
particular, the extremality of the channel is equivalent 
to {|V'n)('0n|} being an operator basis. A constructive 
way of testing may be done using the Bloch representa- 
tion. Recall that to a given normalized complex vector 
IV'n) € we can assign a corresponding generalized real 
Bloch vector 6„ G [l3]. Let cti, . . . ,crrf2_i be or- 

thogonal generators of SU(d), that is, the (Ji are hermi- 
tian, traceless operators obeying tvaiUj = 25ij. Together 
with the identity operator 1 these form an orthogonal ba- 
sis of all linear operators in d dimensions, and we arrive 
at the Bloch representation by defining |V'n)(V'n| —'■ Bn-ff, 
where _B„ = i(|, 6„ ) e K,'' and a = (1, cti, . . . , (Tr2^i). 
For a set of d^ projectors {|V'n)('0n|} forming an operator 
basis, {Bji} is a linear independent set spanning IR'' . 

We thus arrive at the following equivalence (see 
App.ll): 

The channel defined via the dynamical vectors 
, . . . , {ipd^}} is an extremal channel 

^ (15) 
Vol(6i,...,6d2) ^0. 

We are thus able to link the extremality of the phase- 
damping channel to the volume Vol(6i, . . . , 6^2) := 



VI. SYMMETRIES AND 
DECOHERENCE-FREE SUBSPACES 

In qubit systems it may happen that environmental 
influences affect different qubits in the same way. If, for 
instance, the wavelength of a fluctuating field is much 
larger than the separation of the qubits certain qubit 
states accumulate the same random phase and coher- 
ence among such states is preserved. To give an example 
consider a classical fluctuating magnetic field that cou- 
ples identically to all qubits via B{t)J2i'^l =• B{t)T,z- 
In such a case all superpositions of states from an 
eigenspace of will not suffer from decoherence [23| . 
Such decoherence-free subspaces (DFS) can be identified 
in experiments [24^ . 

In the quantum channel formalism the DFS appear 
naturally through symmetry considerations. Assume, 
for simplicity, an N-qubit setup where all qubits are af- 
fected by the environment in the same way. Formally, 
this amounts to the invariance of the channel under per- 
mutations of the qubits. In turn, the set of dynami- 
cal vectors |a„) has to be invariant under qubit per- 
mutations. We conclude that |a„) — |a,„) whenever 
{n\ \n) = {m\ |m). Thus, only N + 1 different dy- 

fN 

namical vectors \bk) occur with a degeneracy of ( ^ 

(the dimension of the corresponding DFS), summing up 
to the total of 2^. 

To give an example, for a two-qubit system with full 
qubit symmetry 1 -H- 2 the most general phase-damping 
channel is made from only three dynamical vectors \ai) — 
\bi) , [02) = las) = I62) , 104) = Ih), such that 



D = 



1 (62I61) (62I61) (63I61) 

(61I62) 1 1 (63I62) 

(61I62) 1 1 (63I62) 

.{hlh) (62I63) (62I63) 1 



(16) 
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and the space {|01) , |10)} is a DFS. These considerations 
can of course be adapted to cases of partial symmetries 
of the environmental influences. 

VII. CONCLUSIONS AND OUTLOOK 

We study phase damping (pure decoherence) from an 
environmental perspective. Any given phase-damping 
channel may be understood in terms of an overlap of 
dynamical vectors |a„) characterizing the channel. For 
a quantum environment these are relative environmen- 
tal states. We investigate how the nature of a phase- 
damping process inflicts with properties of these dynam- 
ical vectors. 

For a single qubit, we infer that any possible phase- 
damping channel is indeed Markovian, that is, a snap- 
shot of some time-continuous Markovian evolution. For 
a single qutrit, we find examples of channels that are 
not Markovian: we give a class of channels we show is 
not infinitely divisible. Remarkably, it turns out that in 
case of Markovian phase damping in arbitrary dimension 
the dynamical vectors may be chosen to be multi-mode 
coherent states. 

For a single qubit a phase-damping channel is of RU 
type. For Hilbert space dimension TV > 4 we discuss a 
physical model of phase-damping dynamics that has no 
RU representation. We find that for a phase-damping 
channel acting on a d? dimensional quantum system, the 
RU property may be linked to a {cP — l)-dimensional 
volume. In a previous article, a link between this (abso- 
lute) volume and the norm distance between the channel 



and the convex hull of unitary transformations was found 

m- 

Our considerations are of relevance for process tomog- 
raphy ^25!] where it is a great challenge to reduce the 
dimension of the parameter space of the process. It is 
clear that any additional assumption about the nature 
of the process (phase damping, RU, Markovian) leads to 
further constraints. Our results allow for a characteriza- 
tion of the channel with a minimal number of parameters 
and should help to speed-up the optimization procedures 
involved [26l |. 
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Appendix A: Diagonal Lindblad Form 

In this appendix we show that the diagonal form of 
the generator C is enough to guarantee the Hamiltonian 
and the Lindblad operators to be diagonal as well. In 
a given basis let the generator of the semigroup 

At = e-^* be diagonal, that is, C{\m){n\) — Zm„|TO)(n|. 

With H = J2rn7i^rnn\m){n\ and Li = Z)mn 

this implies 



z S S - -i(h S -5 h ^ + i V J VZ^'^Z''^(S - V/**^/''^(5 U 

i y k k ) 



Letting m = n, r = s, and n ^ s we see that ||Zmr|P = 
for m 7^ r, so that Imr — Smr^r- Insertion into 
(|Aip then implies h„ir — Smr^r, so that Hamiltonian 
and Lindblad operators are found to be diagonal. In ma- 
trix representation, the generator may thus be written 
as 

Zmn = —ii^m — ^n) + {ln\lm) ~ xdl'mlP + ll'rilP)- 



Itpi) , . . . , \ipd^) define projectors giving an operator basis 
iff the real vectors Bi, . . . , B^p are linearly independent. 



Appendix B: Extremality Criterion 



In order to see the equivalence in Eq. (jl5p we 
have to perform some matrix algebra. The vectors 
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which is the case for [27 
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where Vol(6i , . . . , 6^2 ) denotes the volume of the paral- 
lelogram spanned by the real, (P — 1 dimensional vectors 
{bi,...,bd2}. 
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